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1-^ ' Abstract. - Hydrodynamic equations for an isotropic solution of active polar filaments are 

t~| , derived from a microscopic mean-field model of the forces exchanged between motors and fil- 

Q ' aments. We find that a spatial dependence of the motor stepping rate along the filament is 

O I essential to drive bundle formation. A number of differences arise as compared to hydrodynam- 

ics derived (earlier) from a mesoscopic model where relative filament velocities were obtained 
on the basis of symmetry considerations. Due to the anisotropy of filament diffusion, motors 
^ , are capable of generating net filament motion relative to the solvent. The effect of this new 

U^ ■ term on the stability of the homogeneous state is investigated. 

(N 

o, 

\l , Introduction. - Soft active systems are a new and exciting class of complex fluids to 

which energy is continuously supplied by internal or external sources. Biology provides many 
examples of such systems, including cell membranes and biopolymer solutions driven by chem- 
ical reactions, living cells moving on a substrate, and the cytoskeleton of eukariotic cells [1]. 
The cytoskeleton is a complex network of long filamentary proteins (mainly F-actin and mi- 
'^ ' crotubules) cross-linked by a variety of smaller proteins [2,3]. Among the latter are clusters 
Ph , of motor proteins, such as myosin and kinesin, that use chemical energy from the hydrolysis 
O ' of ATP to "walk" along the filaments, mediating the exchange of forces between them [4-6]. 
The self-organization of motor-filament mixtures has been studied by in vitro experi- 
ments [4-6]. Complex patterns, including asters and vortices, have been observed as a function 
of motor and ATP concentration in a confined quasi- two-dimensional geometry [5,6]. The high 
^ • frequency mechanical response of active filament solutions has also been studied both experi- 
5^ , mentally and theoretically [7,8]. The study of these simplified model systems should lead to a 
better understanding of the formation and stability of more complex structures of biological 
relevance, such as the mitotic spindle formed during cell division [1,9]. 

There have been a number of recent theoretical studies of the collective dynamics of rigid 
active filaments. First and most microscopic, numerical simulations with detailed modelling 
of the filament-motor coupling have been shown to generate patterns similar to those found 
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in experiments [5,6]. A second interesting development has been the proposal of 'mesoscopic' 
mean-field kinetic equations where the effect of motors was incorporated via a motor-induced 
relative velocity of pairs of filaments, with the form of such velocity inferred from general sym- 
metry considerations [10-12]. Finally, hydrodynamic equations have been proposed where the 
mixture is described in terms of a few coarse-grained fields whose dynamics is also inferred 
from symmetry considerations [13-19]. Recently, we established a connection between the 
mesoscopic and hydrodynamic approaches by deriving hydrodynamics via a coarse-graining 
of the kinetic equations [20]. Both the mesoscopic and hydrodynamic approaches share, how- 
ever, an important shortcoming. The rate and strength of the motor-induced force exchange 
among the filaments is controlled by phenomenological parameters whose dependence on mo- 
tor activity is not known. The richness of the phenomena exhibited by the cytoskeleton leads 
naturally to the question of how much of the behavior is specific and how much is generic. To 
answer this question it is important to make the connection between microscopic models and 
'generic' hydrodynamic approaches. 

In this paper we present a first attempt at deriving the motor-mediated interaction between 
filaments from a microscopic description of the forces exchanged between the motors and the 
filaments. Our work establishes the connection between the hydrodynamic equations and the 
microscopic motor dynamics. In particular, it shows explicitly that in mean-field models of the 
type considered here spatial inhomogeneities in the motor stepping rate along the filaments 
are essential to drive bundle formation. 

Kinematics of filament pair. - The filaments are modelled as rigid rods of length / (here 
I should be thought of as the persistence length, rather than the actual filament length) and 
diameter h «l. Each filament is identified by the location r^ of its center of mass and a unit 
vector fii pointing towards the polar end. The mobile crosslinks are formed by small aggregates 
of molecular motors that exert a force on filaments by converting chemical energy from the 
hydrolysis of ATP into mechanical work. Each motor cluster is assumed to be composed of 
two heads, with the i-th head [i = 1,2) attached to filament i at position r^^ = r^ + hiSi, 
with Si the position along the filament relative to the center of mass, —1/2 < Si < +1/2. 
The motor cluster has size Im << I- Hence the attachment points satisfy r^ ~ Tj , or 
^ = r2 — ri ~ sifii — S2'n2. A schematic is shown in Fig. ^ Motor heads are assumed to 
step towards the polar end of filaments at a known speed, u(s), which generally depends on 
the point of attachment. Both filaments and motors move through a solution. We assume 
that the filament dynamics is overdamped and the friction of motors is very small compared 
to that of filaments. Momentum conservation then requires that in the absence of external 
forces and torques, the total force acting on filaments centered at a given position be balanced 
by the frictional force experienced by the filament while moving through the fluid. 

We consider a pair of filaments cross-linked by a single motor cluster. Due to the action 
of the motors, filaments 1 and 2 acquire center-of-mass velocities Vi and V2 and rotational 
velocities wi and 0^2 about the center of mass. We evaluate these velocities in terms of the 
known motor stepping rate, u{s), and of the filaments' orientation. In the absence of external 
forces, any force or torque generated by active crosslinks on one of the filaments of the pair is 
balanced by an equal and opposite force or torque acting on the other filament. This requires 

C,j(ni)uij = -0j(n2)w2j , (1) 

Cri^l ^ -Cr'^2 , (2) 

where Ciji^) = C|pi"-i + Cj-i^ij ~ fi-ifi-j) is the friction tensor of the rod, with C|| and Cj^ the 
longitudinal and transverse friction coefficients, respectively, and C,r is the rotational friction. 
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Fig. 1 - (a)Two filaments connected by an active crosslink. (b)The geometry of the overlap: the 
filaments' centers are separated by ^ = siiii — S2n2. (c) The profile of the motor stepping rate. 



The a-th motor head {a ~ 1,2) at position Tma = Tq + SaH-a has velocity Vma = Vq + 
uo{sa)'a.a + Sai^a X &a, where Sa = uoisa) and ha ~ i-jJa X &a- Motors may also rotate 
relative to the filament at a rate u^ma = '^a + (^ l)"^^^Q(sa)(ni x n2)/|ni x 62], with Oa{sa) 
the angle the a-th motor head makes with the filament to which it is attached. If the motor 
heads are rigidly attached to each other we must have Tmi = i"m2 and w^i = <^m2- Since the 
motors are assumed to be point-like in size on the scales of interest, we neglect motor rotation 
below, i.e., 9a{sa) = 0, which then requires iJi = LO2 = 0. The anisotropy of the rods' friction 
tensor, allows for both relative and net translation of the filaments induced by the action of 
the motors. Using Eqs. ^, we find that the relative velocity v = vi — V2 of the filaments is 



V = uo{s2)n2 - uo{si)hi . 

Equations ^ and Q are readily solved for the filaments velocities, Vi^2 = 
velocity V = (vi + V2)/2 of the center of mass of the pair is given by 

V= A(cr, fii • na) I (1 -2<T)(ni • 62) n2Uo(si) + niMo(s2) 
- 1 - (T - cr(ni • 62)^ niMo(si) + n2Wo(.S2) |, 



±v/2 



(3) 
V. The 



(4) 



where A = (cr/2)[(l - af - a^{hi ■ h2f] , with cr = (C± - C||)/(2C_l). The fact that V ^ 
indicates that motor activity can induce a net motion of the pair relative to the solution. This 
arises from hydrodynamic effects due to the anisotropy of the friction tensor and vanishes 
when C,± ~ C||. Also V vanishes identically for 62 = ±ni, so that V = in one dimension. 

Derivation of hydrodynamics. - The concentration c(r, n, t) of filaments with center of 
mass at r and orientation n at time t satisfies a local conservation law, 

dtc^-v-3-n-j, (5) 

with TZ — nx Vfi. The translational current density, J, and rotational current density, i/, are 

J, = -DuV.c - ^^c V,K. + Jr* , (6) 
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where Dij ~ ksTactCij = Duhifij + D±{5ij — fiifij) is the translational diffusion tensor, 
with £'||,± = kBTactlQ\\,i_ and D^ — ksTact/Cr- The effective active temperature, Tact niay 
be higher than the ambient temperature T due to the motor activity [8]. The potential 14x 
incorporates excluded volume effects which give rise to the nematic-transition in a hard-rod 
solution. It can be written as fcsT times the probability of finding another rod within the 
interaction area of a given rod, [21] 

Fex(ri,ni) ^^ I dv2 I dh2 I I ,5(rf -r2^)c(r2,n2,t)i?(ni,n2) , (8) 



where /^ ... = h'^ ^ i'-in '^^■■- ^'^^ ^ (111:^2) = y^ ^ (^i ' ^2) ■ 

The active current of filaments with center of mass at ri and orientation along fii is 

J^'='(ri,fii)= /dr2 /rffis / / ?9(ni,fi2)m(rf)vi(si,S2,ni,fi2)(5(rf -r2^) 

J J J Si J S2 

xc(ri,fii,t)c(r2,n2,i) , (9) 

where m(r) is the density of motor clusters, evaluated at the point of attachment. Finally, 
Vi (si , S2, fii , n2) is the velocity that filament 1 acquires at the point of attachment of the motor 
cluster due to interaction with filament 2, when the centers of mass of the two filaments are 
separated by ^ = r2 — ri ~ riisi — n2S2, as obtained earlier. Using the (5-function to carry 
out the integration over r2, and assuming a uniform density of motor clusters mg, we obtain 

J'''=*(ri,ni) = ?7io Mn2 / / i?(ni,n2) Vi(si, S2, ni, n2)c(ri, ni,i) 

J J Si J S2 

xc(ri +nisi - n2S2,n2,f) . (10) 

The effect of fluctuations in m(r) will be discussed elsewhere. 

To describe filament dynamics on scales large compared to their length /, we expand the 
concentration of filaments c(ri + ^, 112) near its value at ri. After inserting this expansion in 
Eq. H1U|) . the integrals over si and S2 can be carried out term by term and we find 

Jf'=*(ri,ni) = mo / dfi2 c(ri, fii, t)|('Uii), ,1,^2 + (wiiCj>si,s2Vij 

+ ^{vu^3^k)sus2'^i3'^ik + ■■■jc{ri,n2,t) , (11) 

where {vu^j£,k---)si.s2 — J J ^vii^j£,k--- are moments of the local filament velocity. These 
are expressed in terms of moments of the stepping speed, w'") = J_^,^{ds/l){s/l)"'u{s). If 
u(s) is constant, all odd moments, which control the bundling instability of the homogeneous 
state discussed earlier [11,20], vanish. In mean-field models of the type considered here an 
inhomogeneous velocity profile is crucial to obtain bundling. Spatial variations of u{s) may 
for instance arise from motors slowing down as they approach the polar end of the filament 
due to crowding. This is incorporated here by using the step- like speed profile shown in Fig.^ 
where u{s) is constant along the filament, but vanishes in a small region of extent Im << I 
at the polar end. Using this form, we estimate u^^"^ ~ uq and u^^""*""^^ ~ —uo{lm/l)- All odd 
moments are negative as filaments slow down as they approach the polar end. 

As in Ref. [20] we now assume that on large scales the system dynamics can be described 
in terms of a local filament density p(r, t) and a local filament polarization, p(r, t), defined as 
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the first two moments of the distribution c(r, n, t), 

p(r,t) = / c(r,n,<) , 

p(r,<)- /nc(r,n,t). (12) 

Hydrodynamic equations for these coarse-grained densities can then be obtained by writing an 
exact moment expansion for c(r, n,t) and truncating it at the second moment. The integrals 
over the unit vectors fii and n2 are carried out by approximating the prefactor A of Eq. @ 
and the excluded volume term i? with their mean value over the range of integration. This 
approximation does not affect the structure of the hydrodynamic equations, but only the 
numerical values of the coefficients. 

The nonlinear hydrodynamic equations for the filament density and local orientation of a 
mixture of arbitrarily-shaped filaments and motors in d dimensions will be given elsewhere. 
Here we restrict ourselves to the specific case of long thin rods in d = 2, where D|| = 2D± = D. 
We also linearize the equations around the homogeneous isotropic state, by letting p = pQ+Sp 
and p = (5p. The linearized equations in d = 2 are 

dtSp= IdwHp-^wHp + ^v-p 



dtPi 



192 V 24^/2/ 384 



-A. P^ + Idv'p^ + ^ V. V . p + ^ (l + ^) V,^p 



5Ppol 



s„;2 



384 V 24V2 



1-:T7^^»V'^/^' (14) 



where po = Po'i^o and 



a = -mowou*-^^ ~ moVQUolm/T^ , (15) 

f3 — TTioWo^t — 2moVoUo/'n' . (16) 

The parameter a has the dimensions of a diffusion constant and describes filament bunching or 
bundling, which, in contrast to conventional diffusion, tends to enhance density fluctuations. 
The coefHcient /3 is a velocity and describes the rate at which motor clusters sort or separate 
filaments of opposite polarity. If the motor stepping speed u{s) is constant, independent of the 
position s along the filament, then a = 0. In general, even when a 7^ 0, we expect a << (3. 

The hydrodynamic equations obtained here by using a microscopic model for the motor- 
induced filament velocities differ from those derived in Ref. [20] where these velocities where 
written down on the basis of symmetry in two important ways: 

1. The density equation contains a new term ~ /3V • (pp) that describes filament convection 
along the direction of local alignment. This term vanishes for isotropic objects with 
C_L = Oi ■ It arises for rods because in this case the motors can generate net motion of 
the filaments relative to the solution (V 7^ 0). 

2. The second important difference is that there are no a- type contributions to the active 
currents proportional to gradients of the local polarization. These leads to the absence of 



EUROPHYSICS LETTERS 



ap 




Stationary / ' tt— p — 

^ I ■' / oscillating 





stable 



PP 



Fig. 2 - The phase diagram in the a — (3 plane showing the different regimes calculated from a low k 
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expansion of A+ and the eigen-vector e , ' for po = 0.9. The insets show the unstable modes. 



terms of the form ~ Vi[pjV jPi + ...] in the density equation and of terms ~ V jp[V jPi + ...] 
in the polarization equation. These terms are allowed by symmetry and are indeed 
present in hydrodynamic models written down on the basis of symmetry consideration. 
[11,20] They vanish, however, upon coarse-graining, when averaging over the length of 
the filaments. The absence of such terms does not change the bundling instability of 
density fluctuations, but it does eliminate the possibility of the oricntational instability 
of polarization fluctuations observed in [20] at high filament concentrations. 

The absence of terms proportional to polarization gradients in the active current is not related 
to the fact that our model does not incorporate any motor-induced filament rotation. It arises 
instead from the fact that the motor dynamics is described in terms of the rate u(s) at 
which each motor head steps on each filament. These terms then average out when the force 
exchanged between the two filaments of a pair is averaged along the length of each filament. 

Stability Analysis. - To determine the linear stability of the isotropic, homogeneous state, 
we expand density and polarization fluctuations in Fourier modes, 5p{r) ~ /^ e^^'^pk, p(r) = 
/^, e^'^^'pfc. Defining the longitudinal polarization mode p^ — k ■ pk and V = {dp,pj^), Eqns. 
i)13ll4|l are equivalent to a matrix equation dtV = M. ■ V . The instability of the system to 
density and longitudinal polarization fluctuations is determined the range of parameters where 
for some range of k- vectors, the largest of the eigenvalues A±(k) of the matrix M(fc, po, a, (3) 
has a positive real part. Expressing lengths in units of I and time in units of P / D we define two 
dimensionless parameters a ~ a/D and j3 ~ 1/3 /D [22]. For /3 = the equations are decou- 
pled. In this case there is a stationary instability of density fluctuations for a > ct^ = 3/(2po) 
and < fc < fc°, with A:° ~ (a — 5°)-^/^. A finite, but small value of the sorting rate /3 
tends to stabilize the homogeneous state. This is easily seen from the linearized equations 
as for times 3> D~^ longitudinal polarization fluctuations simply follow density fluctuations, 
and to leading order in the gradients pj^ ~ [Ppo/{2Dr)]ikpk. Eliminating the longitudinal 
polarization in the equation for the density, we find that a finite sorting rate /? effectively 
enhances diffusion, shifting the instability to ap0) = 3/(2po) + (/3^Po(l + l/(4V^))/(24v^). 
The range of wavevectors over which the modes are unstable is again < fc < fc^ (/3) , with 
kaiP) ~ [a — ap{P)Y'^. For larger values of /3 the modes are complex conjugate, describ- 
ing propagating fluctuations. The homogeneous state becomes unstable via an oscillatory 
instability. The location of the instability is still, however, controlled solely by the bundling 
coefficient a. It occurs at the /3-independent value, dm ~ 7.6/po above which Re{X±{k)) has 
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a positive maximum, for an intermediate range of wavenumbers. The actual boundaries of the 
various regions have been evaluated numerically and are shown in Fig. 2. To establish which 
order parameter (density, polarization, or both) goes unstable, we have evaluated the right 
and left eigen-vectors associated with the unstable mode, defined by M • e^{k) = A±(fc)e^(fc) 
and e^{k) ■ M = A±(fc)e'^(fc), respectively. The unstable eigen- vector associated with the 
stationary instability at ap(/3) is dominated by density fluctuations at all length scales. The 
unstable eigen-vector associated with the oscillatory instability at cxm is always a mixed one 
which crosses over from density at small k to longitudinal polarization at large k. 

In summary, the homogeneous state is always driven unstable by a sufRciently large value of 
the bundling constant a. The polarization sorting rate /? tends to stabilize the homogeneous 
state. A large /? also changes the nature of the instability from stationary to oscillatory. 
Finally, an important consequence of the absence of ~ a-type contributions to the active 
currents proportional to polarization gradients is that our model yields no pure polarization 
instability at large length scales. In particular, fluctuations in the transverse part of the 
polarization (~ V x p) that may control the onset of vortex-type structures, are always stable 
to linear order. A numerical solution of the full nonlinear equations is, however, needed to 
establish the precise nature of the inhomogeneous states. 
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